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Abstract 

1^ , The XXZ open spin chain with general integrable boundary conditions is shown to possess a g— deformed ana- 

' logue of the Onsager's algebra as fundamental non-abelian symmetry which ensures the integrability of the model. 
00 This symmetry implies the existence of a finite set of independent mutually commuting nonlocal operators which 
form an abelian subalgebra. The transfer matrix and local conserved quantities, for instance the Hamiltonian, 
are expressed in terms of these nonlocal operators. It follows that Onsager's original approach of the planar Ising 
^ I model can be extended to the XXZ open spin chain. 
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^ ■ 1 Introduction 

Quantum integrable models have always attracted attention, starting from Bethe's work 1 on Heisenberg isotropic 
^ , spin chain, and Onsager's solution of the planar Ising model They possess several important domain of 
' applications which range from (super)string theory to condensed matter systems, and allow to obtain exact and 
nonperturbative expressions for physical quantities (scattering processes, correlation functions,...). Furthermore, 
^ j they are deeply related with nice and rich mathematical structures (like the theory of quantum Lie algebras, 
g— orthogonal polynomials, partial differential and g— difference equations,...). In recent years, much interest has 
been focused on quantum integrable models with boundaries 3, 4 for which several exact results have already 
Ci I been obtained although a large number of questions remain. For instance, the underlying symmetry of these 
models, exact spectrum and correlation functions are yet to be found. 

Among these models, one finds for instance the XXZ anisotropic open spin chain with general integrable 
boundary conditions which has received a lot of attention in the past years [HIIZIEI- Despite of - nevertheless 
important - partial results (some symmetries for a special choice of boundary parameters, spectrum and Bethe 
ansatz equations for special relations between left and right boundary parameters, a connection with the Ising 
model at special values of the anisotropy parameter and special diagonal boundary conditions I^), its underlying 
symmetry for arbitrary boundary parameters clearly needs to be clarified. Indeed, for special diagonal boundary 
conditions at both ends of the spin chain the model is known to enjoy {7^1/2(3^2) symmetry Also, for special 
boundary conditions (left identical non-diagonal and right identical diagonal) it was observed in jllj that some 
nonlocal operators^ act non-trivially on the double-row monodromy matrix (also called Sklyanin operator) and 
Hamiltonian of the open spin chain. Despite of these facts, for non-diagonal boundary parameters neither the 
non-abelian algebra on which these nonlocal operators close (the boundary analogue of [/gi/2 (5/2)), nor its abelian 
subalgebra (generating all integrals of motion and implying the integrability of the system) have been explicitly 

^ baseilha@phys.univ-tours.fr 
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•^Nonlocal in the sense that they act simulateneously on all sites of the chain. Note that the structure of these operators is consistent 
with the one derived in the context of quantum field theory for the boundary sine- Gordon model [T2llT3| . 
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identified. This might be even more surprising, as the term "boundary quantum group algebra" corresponding to 
such algebraic structure is used repeatedely but with no reference to proper defining relations independent of the 
spectral parameter u, derived from the reflection equation^. However, a breakthrough recently emerged [TUfT^] 
from the study of the algebraic structure encoded in the reflection equation for certain (left purely non-diagonal and 
right diagonal) boundary conditions and Ugi/2{sl2) trigonometric i?— matrix. Indeed, in |14| ll5 j it was discovered 
that the boundary analogue of C/^i/2(s^2) is associated with the tridiagonal algebra, first introduced by Terwilliger 
[TB). This algebra is defined through a pair of tridiagonal relations; these are extensions of the g— Serre relations 
and can be viewed as "g— deformed" Dolan-Grady relations due to their close analogy with the relations discovered 
in JJj. In addition, later on it was argued jlH] that this algebra is related with a g— deformed analogue of the 
Onsager's algebra^ constructed in JH|- As briefly mentionned in this remarkable tridiagonal algebraic 

structure survives for the most general (non-diagonal) solutions of the refiection equation. 

In this paper, we show that the homogeneous XXZ open spin-^ chain with the most general integrable 
boundary conditions can be built from a finite set of nonlocal operators which generate a g— deformed analogue of 
the Onsager's algebra and ensure the integrability of the system. The symmetry of the model is identified. Then, 
our results open the possibility to apply Onsager's original approach [2j to integrable open spin chains. 

The paper is organized as follows. In the next Section, it is shown that the entries of the double-row monodromy 
matrix (Sklyanin operator) can be written in terms of a finite set of nonlocal operators acting on finite dimensional 
tensor product representations and generating a g— deformed analogue of the Onsager's algebra. Then, we study 
in details the abelian symmetry of the model: the Hamiltonian is shown to commute with a finite set of nonlocal 
operators which generate an abelian subalgebra of the g— deformed Onsager's algebra. In the last Section, 

based on the knowledge of the symmetry we argue that the Onsager's approach can be applied to the XXZ 
open spin chain. Indeed, we explain that the spectral problem for the transfer matrix, and in particular the 
Hamiltonian, reduces to the one for the nonlocal operators ^2k+i' fc G 0, iV — 1. 



2 The XXZ open spin chain with integrable boundary conditions 

The XXZ open spin chain with general integrable boundary conditions has recently attracted much attention. 
This started from the work of Sklyanin 0], who proposed a systematic procedure to find boundary conditions 
for open spin chain compatible with integrability. Since this work, the objective remains the same: derive exact 
results such as the energy spectrum of the open spin chain for arbitrary integrable boundary conditions and find 
exact correlation functions, extending the known results for the infinite volume limit (see for instance 119! and 
references therein). In particular, in the infinite volume limit the identification of 11^1/2(312) as the underlying 
symmetry of the spin chain played a crucial role in deriving such examples of exact results. Similarly, identifying 
the underlying symmetry of the XXZ open spin— ^ chain for general integrable boundary conditions is a very 
important problem, which is completely solved in this Section. In the following, we consider the Hamiltonian 

w-i 



tixXZ = 2^ y(^l(^l +(^2^2 +^'^3f^3 

fe=l 



(,W+,(0)) V 2 " (gl/2-g-l/2 

(gi/2 _^-l/2) /(k-K*) 



■^^2{q''^ + q-'/'){n+a^ + n^a^)) , (1) 



(k -I- K*) 



■'By analogy, the Uq{sl2) algebra is associated with a set of fundamental defining relations independent of u, although derived 
from the ^^RTT^^ relations (Yang-Baxter algebra). 

^For 9 = 1 and special value of p, the algebra introduced in 1181 can be related with the Onsager's algebra which played a crucial 
role in the original solution of the planar Ising model (see for details). 



2 



where cri.2,3 and a± ~ (cti ±i<T2)/2 are usual Pauli matrices. Here, A — {q ' +q ' )/2 denotes the anisotropy 
parameter. Boundary parameters associated with the most general integrable boundary conditions are such that 



(o)^t - 



k+ = {k_y = -{q'r^ - g-i/2)e'V2 



)t = -f^i/s 



-1/2 



)) 



(left) , 
(right) . 



(2) 



In order exhibit independent parameters, the r.h.s parametrization can be replaced in ^ in which case one imme- 
diatly identifies the Hamiltonian as defined in 7.. In total, one has six boundary parameters cqq, cqi, cqo, cqi, 9, 9. 
Considering a gauge transformation, note that one parameter might be removed in which case only five arbitrary 
parameters remain. For symmetry reasons, we however keep the boundary parametrization as defined above. 
Note that restricting the boundary parameters to special values, one obtains the cases considered in [HI5l[TT|. 

A characteristic feature of quantum integrable systems is the existence of an infinite number of mutually 
commuting integrals of motion. For quantum integrable systems on the lattice with non-periodic boundary 
conditions, Sklyanin ^ proposed a rather general method to derive in a systematic way this family of independent 
mutually commuting quantities. Similarly to the periodic case, a transfer matrix which provides the generating 
functional for all integrals of motion can be constructed. For the homogeneous XXZ open spin— ^ chain with N 
sites, this transfer matrix takes the form 0] 

(-1)^ 



txxziu) = 



-l\N 



trn 



(w-^ + u ^ — q — q 

where tro denotes the trace over the two-dimensional auxiliary space, the Sklyanin operator 

K''^\u) = %(uuAr) • • • Li{uvi)K^{u)Li{uVj;^) ■ ■ ■ %(uw^^) 
defines ® a family of solution to the reflection equation 

R{u/v) (if(^)(u)® J) R{uv) (lI^K^'^^v)) = {I®K'-^\v)) R{uv) (K'-'^^u) ® I) R{u/v) 



(3) 



(4) 



(5) 



for arbitrary parameters vi,...,vn e C and tt's) induces a spin— ^ representation for the quantum space located 
at each site of the spin chain. As shown in the integrability relies on the fact that given a i?— matrix the 
Lax operator L{u) satisfies the Yang-Baxter algebra, K-{u) in Q is itself a, solution of |SJ and K^{u) solves 
the "dual" reflection equation^. In the context of integrable open spin chains, the element K^{u) and K^{u), 
respectively, characterize the class of left (resp. right) integrable boundary conditions. For the XXZ open spin 
chain with transfer matrix Q, given the trigonometric i?— matrix associated with ?7 1/2(5/2) 



/ uq'/^^u-^q-'/^ 



R{u) 



7I/2 



7-1/2 



gl/2_g-l/2 
U — U^^ 





\ 



uq^/^ — u^^q^^/^ J 



(6) 



the most general elements K±{u) with c— number entries take the form 



K+{u) = 



-1 (0) 
u e_ 



u~^)/{q'l^ 



7-1/2A 



fc_(M2-u-2)/(^l/2_g-l/2) 
vl/2, 



(0) 



K^{qu — q 



q ' UK 
2 „-l„,-2 



q-^^^u-^K* 



u -){q 



1/2 



7-1/2^ 



ue 
Hi+{qv? 

q 



lu-2)(gl/2 



1-1/2^ 



(7) 
(8) 



with |(5J), respectively. For further convenience, note that we have chosen a different normalization in the off- 
diagonal boundary terms between and 

^Note that K^^\u) acts on a tensor product representation of "quantum" spaces. Standard matrix multiplication over the 
auxiliary space is used to multiply the L— operators. 



''The "dual" reflection equation follows from j^J by changing u 



-1/2, 



^/^i) ^ and K(u) in its transpose. 
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With the definitions above, differentiating the transfer matrix Q n— times with respect to A = ln(?i) at A = 
(i.e. u — 1), mutually commuting local conserved quantities can be derived. Recall that 

Rokiu ^ I) ^ iq'/^ - q-'/^)rok , 

where Vki denotes the permutation operator with properties 

Vki = Vik , Vli = 1 , PkiFi„, = FkmVki m^kj 

for an arbitrary operator Fi^i acting on sites k and I. Then, for n = 1 one derives in a straightforward manner 

^ ,/• i/2_ -i/2\ 2N \ 2 

where the Hamiltonian given by ^ is one of the local conserved quantities. More generally, higher mutually 
commuting local conserved quantities, say Hn with Hi = Hxxz can be derived similarly. 

2.1 Preliminaries: structure of the Sklyanin operator 

The Sklyanin operator, as defined in Q is one of the fundamental element in the transfer matrix of the XXZ 
open spin chain associated with the Hamiltonian JQ). For further analysis, let us write it in the form 

^ w- cW(u) pW(7.) ) ■ 

Remarkably ^H], it is possible to systematically disentangle for arbitrary integer N the algebraic part (acting on 
tensor product representations of N quantum spaces) from the spectral parameter (i.e. dependent on "u") part. 
Indeed, using the procedure presented in ,15 , for Q, ijHJl and arbitrary iV after some rather involved calculations 
(see ^1 for details) the Sklyanin operator H1U|I is such that 

N-l N-1 

^pi-)(«)wLl^-.-V^/^E 

k=0 k=Q 

k=0 fe=0 
,2-„-2Wfc,fc ^ol/2,;2^„-l/2„-2^, T ^1 

k=Q 
N-1 



A^'^Hu) = lief) + u-hi""^ + iu' - 7i-2)(^qV2 ^ p[1\u)\N'^_\} - u-'q~'/' ^ p['l\u)\N^^^\] 

k=0 k=Q 
N-1 N-1 

fc=0 fe=0 

B<".(„, . ^^^( ^-^^-',';;;!;_c;'""' ^rv) + Y./-"M^i^'. * -r- ) 



Compared to although apparently complicated this expression for K'-^^u) possesses rather nice features. It 
contains non-algebraic (c— number valued) terms such as the Laurent polynomials P^_^\u), e^^\ '^o^'* separated 
from "u"- independent operators W_fc, W^+i, G^+i, Gk+i- Namely, one finds 

• Non-algebraic part: The "initial" solution K-{u) of the reflection equation being slightly different 
from the choice in | 18| (here, diagonal boundary parameters e^'' are considered), the final expressions for the 
Laurent polynomials become 



N-1 



P-Vi-^)- („l/2 , „-l/2) E „l/2 , Vl/2 ) (12) 



2 V n—k 



(qi/2 + g-i/2) ^ V gi/2+q-i/2 



N 

^ (_i)A'-n(^l/2 ^ g-l/2y.+ l ^ . . . 

fel<...<fcjV-Tl-l = l 
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where we have defined 



The remaining constant terms in are 



^ ^ /c— 1 

for arbitrary values of A^. In all expressions above, 

^0.) ^^,.+1/2^^-^.-1/2 



N 



denotes the Casimir operator eigenvalue associated with the spin— representation of [7^1/2(5/2) for fc = 1, N. 
Note that setting e^"* — 0, jk = j, A;+ = 1/co, fc_ = 1 and Vk = w in above expressions, one recovers the results of 

CHI 

• Algebraic part: For general values^ of N, the elements W_fc, Wi:+i, Gfe+i, G^+i (4iV in total) act on 
iV— tensor product representations of Uqi/2{sl2), and depend solely on the A^— parameters Vk and spin— for 
fc = 1, A^. Following an analysis similar to 18 for (0, ^ one eventually gets the family of nonlocal operators 
which act as: 

-I. (^1/2 _^ g-l/2) k (-^1/2 _|_ g-1/2) fc+1 (^1/2 _|_ ^-1/2)2 W 

/c+l (^1/2 _^ ^-1/2) A.+1 (^1/2 _^ g-1/2) (^1/2 _^ g-l/2)2 k 

^fc+l ~ A:+((jl/2 +q-l/2) ^ (^1/2 _^g-l/2)l."jv9 +t^Ar9 j » + « (a^^+l 

+ (g - q-') (k^v^q-'/^S^q^^/^ ® {W^-'^ W^^-^)) + k^y-^V^' S^q'^'^^' ® (W^^^^^ - "^^l+l^)) 

(gl/2+g-l/2)2 fe ' 

'^fe+l - fc_(ql/2 +^-1/2) (^1/2+^-1/2)'^^^^'? +^'Ar9 j » ^fe + ^ <^fe+l 

+{q Q-') {Kv-^\"'s+q^^'^ ® (wL^-^^ - wf -^)) + fc+«^g-i/45+g-3/2 ^ (vvi^-i) _ wLl;f )) 
, « + .... 

(gl/2 +q-l/2)2 '^fc ' y"^^) 



*To understand the general case A'', it is instructive to look at the special cases TV = 1 and N = 2 first 1151 1181 . For Af = 1 in 
in . it is easy to show that Wo,Wi are simply given by linear combinations of 1/^1/2(^/2) generators SiQ^^^/z g±s3 with coefficients 

depending on the parameter v\, whereas Gi, Gi are linear combinations of S^, q^'^ and S±q^''^/'^ (for non-vanishing values of e|^'). 
For N = 2, the elements Wq, W_i , Wi , W2 and Gi, G2, Gi, G2 are linear combination of tensor product of U^i/2{sl2) generators, with 
coefficients depending on vi,V2. 
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where, for the special case k = we identify ^ 

Wj. |fe=o = U , W_j._^Jfc=o = U , G^. |fc=o - ^fe |fc=0 = (gl/2 _ g-l/2) ^ ■ '■^^^ 

Here, the identity operator j(^) = J (g) • • • ® J acting on the A^— quantum spaces tensor product has been 
introduced. In addition, one has the "initial" c— number conditions 

w(")ee6V"', w'^^EEeL"' and gJ") = gS"^ ^ 6^ 6^"^ (qi/2 _ ^-i/2) . (^g) 

For the sake of complctncss, let us mention some important properties of the elements H14|l given above, 
which can be deduced inserting K^-^\u) in the reflection equation (jSJ. Expanding in the spectral parameter 
u, one obtains a pair of relations linear in terms of the elements W_fc, W^+i which ensures the closure of the 
algebra. For = 1, they correspond to the Askey- Wilson relations |14[ I15|. trilinear in terms of Wo,Wi or, 
equivalently, linear in terms of Wq, Wi, W_i, W2 in which case W_i, W2 are defined as nonlinear combinations of 
Wq, Wi. For N = 2, one obtains relations of fifth-order in Wo, Wi, or linear in terms of the fundamental elements 
W_fe, Wfe+i, /c = 0, 1, 2 having identified W_2, W3 with certain nonlinear combinations of Wq, Wi ^H). In general, 
given N a pair of relations of order 27V + 1 in Wq, Wi exists, which can be written as a pair of linear recurrence 
relation of length A^-|- 1 in the fundamental elements \N^k,^k+i, fc = 0, 1, ...,7V where W_Ar, W^r+i are nonlinear 
combinations of lower elements W_fc, Wfc+i with A; = 0, 1, — 1. In the present work, for the choice lO, lO 
straightforward calculations similar to ^Hl yield to 



N 



k+k- 



k+k- 



k=l 

N 

(N)^^r(N) . > 



W(N) ^ ^ Ci^^^ + J(^) = . (17) 



fc=l 

Also, linear relations for the elements G^+i, Gfe+i can be derived using ifTH) in lfT7|l . They read 



N 



+ fc=l 

-^^L^^f^j-W^ + E ^'-%A1\ - . (18) 
+ fc=i 

Having the Sklyanin operator in the form pU|) with let us apply above results to the XXZ open spin 

chain. Plugging H10|l with Ijllll in Q, the Lax operators in Q which characterize the contribution of the "bulk" 
site— fc spin degrees of freedom in Q reduce to i?— matrices of the form ^ and denoted above RQk{u). As a 
consequence, the expansion of the transfer matrix at u = 1 automatically leads to local conserved quantities 
0]. Setting jk = 1/2 and Vk = I for all k = 1, ...,iV (homogeneous case), the Sklyanin operator associated with 
the XXZ open spin chain which enters in JSJ leads to consider the nonlocal operators defined by 

wiT' - («f.i-^^^)[wLT^]k^i , >v(f,'^(®-,.(^))[w|^,']k.,, 

Gi:"! - «.i-^^))[G|^i^]k.^i , C^-«^i-^^^)[C]k.^i- (19) 
for I S 0;...,A'^ — 1. Using the two-dimensional representation 

^(1/2) [^^] = and ^(1/2) [gg] ^ , (20) 

^Although the notation is ambiguous, the reader must keep in mind that W^'^'I^j—q 7^ W^'|;j_q ,W^'^-^ |j,_q 7^ W^^'^'j^ |j._o for 
any A'^. 

^'^More general relations can be obtained as a consequence of these fundamental ones. For details, we refer the reader to 1181 . 
Appendix B. For (J = 1, note that the relations 1171 . 1181 as well as their generalizations 1181 reduce to the finite recurrence relations 
for the generators A^, Gi in Onsager's algebra |21|. 
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their explicit expressions can be easily obtained from H14|l . For instance, the simplest nonlocal operators take a 
rather simple form: 

W^"^' = (fc+a+ + fc_a_)® 1(^-1) +g-'^^/^®>V^'^-'^ (21) 

with "initial" conditions 

^ - k+a+ + k^a^ + e^^\--'/^ . 

Consequently, the Sklyanin operator I^J associated with the homogenous XXZ open spin chain takes the form 
pU|) with where the non-abelian algebraic structure - originally hidden in the reflection equation Q - is now 
encoded in the nonlocal operators p9|) . For special boundary conditions, note that the nonlocal operators (|21|l 
coincide with the ones obtained in jllj [l^ . 

2.2 Non-abelian symmetry: q'— deformed Onsager's and tridiagonal algebras 

The integrability of the XXZ open spin chain is a consequence of the reflection equation |SJl or, more precisely 
its underlying non-abelian symmetry. Indeed, as shown in 4 transfer matrices of the form Q are mutually 
commuting provided K^^\u) solves (O. As the algebraic content of the Sklyanin operator K^'^\u) is now solely 
encoded in the finite set of nonlocal operators H19|l . our next task is to identify the algebra on which they close. 
As a consequence, this non-abelian algebra ensures the integrability of the XXZ open spin chain. 

To exhibit non-trivial commutation relations among the nonlocal operators H19|) , one can consider the asymp- 
totic expansion around u — > oo of the reflection equation ^ for an arbitrary spectral parameter v. In this limit 
one easily gets from Q and (|10|l with Hll|) the simple intcrtwinning relations ^1 

Using the explicit expressions H14(l together with H20() . explicitely one finds 
(7r(V2)xzdW)[Wr+^)] ( 9^Wr^ 

(^(V2)^,rfW)rw(^+i)i ^ / g-^w(^^ k+v-^ \ 

It is then easy to identify the symmetry behind the reflection equation |14j . As an intermediate step, the following 
algebraic relations among the entries of are immediatly obtained from l|22ll : 

[\N^^\A^^\v)] =g-l/2i;-l(/fc_i3W(i;)-fc+CW(«)) , 
[W[,^\PW(1;)] = -q'/\{k^B^^\v) - k+C^''\v)) , 

^ g(iV) („)] ^ ^ („_4(W) _ y-^V^N) ^ 

[W(^\CW(.)]^_, ^ -fc_(.^W(.) - .-ipW(.)) , (23) 

where the commutator [AT, F], — q^/^XY — q^^/^YX has been introduced. Similar relations are obtained 
for \n[^\ provided one substitutes q — > q^^,v v^^ in 123|) . Then, plugging the explicit form of the ele- 
ments A^^''{v), B^-^\v), C^-^\v), V^^^v) in these relations, g— deformed commutation relations for the nonlocal 

^^Setting V = exp(^A), g^/^ = exp{i^) and choosing the special case = fc_ = 1 one recognizes the relations proposed in 1111 . 
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operators ((T^ are easily obtained. It yields to 



;(N) ^{Nh _ r (N) (N). _ (gfc+i - Gk+l) 



[wr,>viT]=0, [>Vr,w(^}]=0 (24) 

and 

for k,l £ 0, N — 1, where 

p^iq'/^ + q-'/^k+k^ . (25) 

Comparing if^ . with (eq. (4)) we find that the operators Vvi^\ ^fc+i' Gk+i^ ^i+i induce a representa- 
tion of the g— deformed analogue of the Onsager's algebra. We now explain how Wq^\ w[^^ give a representation 
of the tridiagonal algebra^^ recently introduced by Terwilliger By Definition 3.9] the tridiagonal algebra 
T is the unital associative algebra generated by two symbols A, B subject to the relations 

[A [A [AS] J,-.] = P[A,B] , (26) 

[bAbAB.A]^-^] = P[B.A] . (27) 

The equations H2t)|) . (|27|l are called the tridiagonal relations or the q-deformed Dolan Grady relations. Based on 
jl4l [TH] and previous results one has 

[wr,[wr,[ww,wr]j,_,] = (28) 

Comparing (|26() - (|28|l we find Wq^\w\^^ induce a representation of T. We remark that for iV = 1,2 each of 



a 



^-fe ' ^fc+i' ^fe+i' ^fc+i ^ polynomial in Wq , W} as shown in JS]. From these comments, we see 
close relationship between T and the algebra (|24|l . The exact nature of this relationship is unclear and needs 
further analysis. Indeed, it might be possible to construct a different basis of nonlocal operators, say Ak, Qi, Gh 
which g— commutator relations would be exactly in one-to-one correspondence with the defining relations of the 
Onsager's algebra. Although finding such basis is an interesting problem, we leave it for future work. 

Finally, let us mention that in the context of quantum integrable systems with boundaries the algebraic 
structure H28|) discovered in |14l I15| provides the first explicit example of boundary quantum algebra. 

2.3 The symmetry of the XXZ open spin chain 

Due to the simplified expression pU|) with (|ll|l for the Sklyanin operator it is now possible to expand the transfer 
matrix in terms of linear combinations of the nonlocal operators p9|l . After some straightforward calculations, 
the transfer matrix ||2J) becomes 



txxz {u) = Yl -^2^+1 (u) + -^0 {u) (29) 

with 



k=a 



^It is a coideal subalgebra of U^i/2(sl2) |14II15| . For higher rank quantum Lie algebras, similar relations can also be derived. 
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In (|29|l . the Laurent polynomials are given by 

[q + q ^ — — u 

= {q^q-^-u^^u'^r V^ + ^H ^0 (") + "0 

((qV2y2 ^ ^-1/2^-2)^ ^ (gl/2 ^ q-l/2)^*)gW 

(q + - u2 _ u~2)Af 

((ql/2y2 ^ ^-1/2^-2)^* ^ (^1/2 ^ ^-l/2)^)gW 
^ (q + g-1 - U2 _ u-2)Ar 

with the identifications 

^^-fe (^) = -P-lfc (^)l{j,=i;«,=i v;} and t^^, ' = ^ ' |{j,= i ;„,=i v;} • 

As mentionned above, the reflection equation implies the commutativity of the transfer matrices of the 
form ||2Jl which is nothing but a manifestation of the abelian symmetry of the XXZ open spin chain. Indeed, using 
the properties of the reflection equation and the Yang-Baxter algebra for arbitrary spectral parameters u, v it is 
not difficult to show that @] 

[txxz{u),txxz{v)\^0 . (31) 

As an immediate consequence, the nonlocal operators X^^ii ^ = 0, N — 1 as defined in %\Vi\ generate an abelian 
subalgebra of the g— deformed analogue of the Onsager's algebra: replacing the expansion in we deduce 
that all nonlocal operators H3U|I are in involution i.e. 

[likli^li^+i] = for aU k,lcQ, N - 1 (32) 

for generic values of the left and right boundary parameters. As a consistency check, note that it is possible 
to derive the same result using solely the defining relations (|24|l . Also, for q = 1 and k± — the hierarchy of 
nonlocal operators 12^1' k = 0, .... N — I coincides with the hierarchy discovered by Onsager in the Ising model 
12 or, equivalently |2(J[ I21| . the one discovered by Dolan and Grady in |17j. 

Finally, combining the expansion of the transfer matrix in terms of local conserved quantities together with the 
expansion in terms of nonlocal ones, we deduce the abelian symmetries of the model. Using ^ in H31() together 
with the expansion ()29|l for u v we end up with 

[Hxxz,4kli] = for all fc e 0, ...,N- 1 . (33) 

Obviously, it is possible to restrict the boundary parameters to special values. For instance, for k* = and 
K± = 0, the Hamiltonian commutes with all nonlocal operators w'^\ k = 0, ...,N — 1. In particular, the case 
fc = was identified in fIT}. 

3 Comments 

We have shown that the integrability of the XXZ open spin chain with general integrable boundary conditions 
relies on the existence of nonlocal operators of the form H19|l generating a deformed analogue of the Onsager's 
Lie algebra with defining relations (|24() . The transfer matrix of the model has been expanded in terms of a 

^''For the homogenous case and spin— ^ representation at all sites, one simplifies 

2(g + g-l) (;v-l)! / 2{g + q-l) TV e^°\q^/^ - q-^/^)2\ 



E 



ail ■ ■ • 



(gi/2 + g-i/2); (fc)!(Ar _ fc- 1)! V(gi/2 +q-i/2) TV _ A; ' fc+fe_(gi/2 + g-i/2) 
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hierarchy of mutuaUy commuting quantities (|30|l which generate the abehan symmetry H33|) of the modeL In 
particular, these quantities - which generahze the ones of Onsager Pj or Dolan-Grady 17 - are expressed as 
hnear combinations of the nonlocal operators (|19|l with coefficients depending on the boundary parameters. 

Clearly, the results presented here open the possibility to study the XXZ open spin chain - and more generally 
quantum integrable models with boundaries - from a new point of view, based on Onsager's approach to integrable 
lattice systems. To motivate the proposal, let us first recall the "Onsager's approach" . In the pioneering work of 
Onsager [2] for the two-dimensional classical Ising model in zero magnetic field, the introduction of the so-called 
Onsager's Lie algebra with generators A^, Gi and defining relations 

[Ak.Ai] = AGk-i , [Gi^Ak] = 2Ai+k - 2A^i+k , [G,-, Gi] = (34) 

for any integers k, I played a crucial role in the analysis. Actually, Onsager was able to show that the transfer 
matrix of the Ising model can be written in terms of the fundamental generators of H34() . For an arbitrary 
parameter A, it takes the simple form^'' (see \'21\ ) 

<7sms(A) = exp(AaoAo) exp(AaiAi) exp((l - A)aoAo) , (35) 

where ao,ai are parameters depending on the temperature of the system and the coupling constant in each 
direction of the two-dimensional lattice. For a rectangular lattice of size L x N (with site index "i" in one 
direction) , the operators can be represented by matrices of dimension 2^ which read |21 121| : 

N-l N 

j=i i=i 

Here the choice of sign in the definition of Aa corresponds to periodic/antiperiodic (-1-/—) boundary conditions 
in one direction of the lattice. For any boundary conditions, it is known that Aq,Ai satisfy the Dolan-Grady 
relations ^7] 

[Ao,[Ao,[Ao,Ai]]]= po[Ao,Ai] and [Au[Au[Ai, Ao]]] = po[Au Aq] , (37) 

where po — 16. Although these relations clearly follow from H34|) . it is important to mention that these relations 
are actually sufficient to generate the Onsager's algebra |20[I21[[23] and imply the integrability of the Ising model 
with transfer matrix H35|l . For generic values of ao, ai, using the Dolan-Grady relations l|37|) in the expansion of 
H35|l it is possible to show [2| that the quantum Hamiltonian 

, ,„ , . , \ ^ sinh(2 — 4A)ao ^ sinh(2Aao) sinh(2(l — A)ao) , , , , 

H = coth 2ai + coth 2aoMi + . [ Gi ^ . \ —{Ai - A_i 38 

smh(2ao) smh(2ao) 

satisfies [t ising, H] ~ 0. More generally, any local conserved quantities can be written as linear combination of the 
fundamental operators A^, Gi which obey H34() . Based on the properties of the non-abelian Onsager's algebra 1)3411 . 
Onsager found the eigenvalues of the transfer matrix in the principal direction which lead to the determination of 
the free energy. Later on, this approach was applied to a few other examples of quantum integrable spin chains 
which enjoy the Onsager's symmetry algebra H34I) - Among them, one can cite for instance the XY model [17| . 
the superintegrable chiral Potts model [221 121[|^I24| and generalizations For each of these spin chains, the 
Hamiltonian and more generally any higher local conserved quantity are simply expressed in terms of operators 
which form the abelian subalgebra of Onsager's algebra. Remarkably, this finite set of linearly independent 
operators coincides with 1)30(1 provided one identifies g = 1 in ()24|l and fix the boundary parameters k,k*,k± to 
suitable values. For all these models, the spectral problem for the Hamiltonian or higher conserved quantities 
reduces to the one associated with the family ()30|l at q — 1, which solution is well-known [2 1^ 1221 ■ 

By analogy, it is thus natural from previous analysis to consider the XXZ open spin— ^ chain with general 
integrable boundary conditions from that point of view. The transfer matrix of this model Q being simply 

^*In original Onsager's paper |2] the parameter A = 1/2 or A = 1. 
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expressed in the form H29|) i.e. in terms of the generators (|30|l of the deformed Onsager's abehan subalgebra 
H32|) . any local conserved quantity can be written in terms of H3U|) differentiating txxz{u = exp(A)) given by (|29|l 
with respect to A at A = 0. For instance, using Q one immediately gets 



Hxxz 



(,l/2_g-l/2) 



) 




(39) 



For small values of iV (TV = 1 m US], = 2, 3), it is not difficult to solve numerically the spectral problem 
for l2k+i ^'-'^ fc = 0, ...,N — 1 which exhibits some interesting features depending on the boundary parameters. 
For instance, for non-vanishing values of the boundary parameters the spectrum of the simplest quantity is 
non-degenerate. However, degeneracies appear for k* — k± — (or equivalently k = k± =0). For higher values 
of N important simplifications in the analysis occur for special relations |51 [7] between left and right boundary 
parameters. In this case, the eigenstates of (|3U|) can be factorized in agreement with the results of Details 
will be considered elsewhere. 
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